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<^ \ Abstract 

! The hadronic screening mass at high temperature (T) in QCD 4 is examined on 

the basis of the QCD sum rules in (2+1) dimensions. Due to the magnetic gluon con- 
densate at high T which is expected to be nonvanishing, the screening mass deviates 
from the asymptotic value 2ttT. Also, the screening mass in the vector (pseudo- vector) 

00 ' channel turns out to be heavier than that in the scalar (pseudo-scalar) channel. 
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Quantum Chromodynamics (QCD) exhibits a transition from the hadronic 
phase to the quark-gluon plasma at finite temperature (T) and chemical potential. 
In fact, numerical simulations on the lattice show a rapid growth of the internal en- 
ergy near T c ~ 150 MeV [1]]. Originally, it was expected that the system above T c 
behaves as a weakly interacting gas of quarks and gluons. However, there have arisen 
several indications against the simple picture at high T, which include the infrared 
breakdown of perturbation theory 0, area law behavior of the space- like Wilson loop 
@, deviation from the Stefan-Boltzmann law for pressure [§J, non- vanishing gluon 
condensates ||, and non-perturbative effects in hadronic correlations || 0]. 

The purpose of this paper is to examine the relation between two non-perturbative 
effects at high T, i.e., the magnetic gluon-condensate and the hadronic screening 
masses. The latter is related to the meson masses in the QCD 3 + adjoint Higgs sys- 
tem which is a high T effective theory of QCD4 (see || for the recent developments). 
We will use the QCD sum rules in (2+1) dimensions to study the spectra of the 
effective theory. Relation of our approach to the previous works |10, [111 will be also 
discussed briefly. 

The hadronic screening mass /i(T) is defined by the space-like hadronic corre- 
lation at finite T, 

n^(J(0,x)J(0,y))~e-^ T )l x -yl (|x-y|-oo) , (1) 

where (•) denotes thermal average and J(r, x) is a color-singlet hadronic operator. 
In the following, we will focus on mesonic correlations and take J = q^oq for the 
vector channel and J = qq for the scalar channel for simplicity. E3 If the system is 
nearly a free quark-gluon gas at high T, /i(T) = 2nT (nT being the lowest Matsubara 
frequency for fermions) should be satisfied in both channels . However, numerical 
simulations on the lattice show a sizable difference between /i(T) in the vector channel 
and that in the scalar channel, which indicates non- negligible correlation above T c . 
This phenomenon could be explained by (i) perturbative gluon exchange, in particular 
short-range spin-spin interaction giving a difference between vector channel and scalar 
channel, and/or (ii) non-perturbative gluonic contribution surviving even above T c . 
We will focus our attention on (ii) in this paper and analyze the hadronic correlations 
under the non-perturbative background fields using the QCD sum rules. Since we 
will start with the high T effective theory ||, our analysis will be applicable only at 
relatively high T. 



^dn the chirally symmetric phase above T c , /^(T) takes the same value in the same chiral multi- 
plet. Also, we will not consider the disconnected diagrams in (|l|) for simplicity. The lattice QCD 
simulations currently available also make this approximation [Q . 
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The fermionic part of the functional integral at finite T under the background 
gauge field reads 

Z(A,A A ) = J[dqdq]exp[-jyrjd 3 x qi^D^q] , (2) 

where 7 matrices satisfy {7^,7^} = —25^ and the color indices are suppressed. All 
the non-perturbative effects are hidden in the gauge sector of the functional inte- 
gral which is not shown here. At high T, the relevant gauge degree of freedom is 
the "static" background A M (0, x) f|. Let us first substitute the following Fourier 
decomposition into eq.(fj): q(r, x) = T^„e ,u " T g(w n ,x) with u n = 2ir{n + \)T. By 
absorbing extra 7^ into the definition of q and rescale the quark fields as q — > VTq, 
one immediately arrives at Z written as a sum of 3-dimensional quarks with "quark 
masses" uj n moving in the background field = (A(x),0(x)): 

Z(A, <P)= f [dqdq] exp[- £ g„(x)(irf Df -u n + #)g n (x)] , (3) 

where we have defined new 7 matrices Tf = —i'Jo'yf satisfying {Fi ,Tj E } = —25^-. 
For later convenience, we have not rescaled the background gauge fields, thus g in 
eq.(HD is the dimensionless coupling constant in QCD 4 . Since A^ is static (r inde- 
pendent), there is no transition among quarks with different n. Thus, at high T, we 
can safely take n = sector to compute the space-like correlation (|l|): Contributions 
from ti^O sectors are exponentially suppressed relative to the n = sector. In the 
following, we will omit the suffix n assuming that n = sector dominates. 

By taking into account the redefinitions of the quark field we have made, eq.([l]) 
in the vector and scalar channel can be written as 

n 5 oc ^ 7 fg(x)gz 7 fg(0))-e- m S l x l, (4) 
ITy oc (qiq(x)qiq(0)) -> e~"V |x| . (5) 

The vector (scalar) correlation in QCD4 at finite T is reduced to a scalar (skew 
scalar) correlation in the 3-dimensional Euclidean space at zero T. Our 4-component 
representation of the Dirac field in 3-dimensions contains two-layered 2-component 
spinors in 3-dimensions. This is why we have two kinds of scalars: qq and q^fq. 
Similar 4-component spinor is adopted in to study the chiral symmetry breaking 
in three dimensional QED. 

/i(T) in ([!]) is identified with the masses m sv defined in (HH). A possible way 
to evaluate m sv is to make an analytic continuation of Hs,v i n t° (2+l)-dimensional 
Minkowski space and look for the lowest energy poles. This is analogous to the "funny 
space" trick adopted in ref.[[L_t], [□]] where the temporal and one space direction is 
interchanged in QCD3. 
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An major assumption in our procedure is that the effective theory is in the 
confining phase and Hs,v has always isolated hadronic-poles in the Minkowski space. 
This assumption seems to be supported by the recent lattice QCD simulations where 
the effective theory at high T is shown to be in the confining phase and not in the 



Higgs phase, namely ((f)) = |13j . Thus, in the following, we will set = in (y) and 
keep only A(x) which is essential for confinement in (2+l)-dimensions. 

To extract the pole positions of Hs,v in the Minkowski space, we are going to 



use QCD sum rules [£L4[] in (2+l)-dimensions which is suitable to relate the spectral 
parameters and the background fields (condensates). Here, (G^) is the relevant oper- 
ator with the lowest dimension and gauge invariance in our analysis. This is nothing 
but the magnetic condensate (B 2 ) in 4 dimension. Lattice simulations of (B 2 ) in 
QCD4 at T > T c show that it is non vanishing at least for T c < T < (2 — 3)T C H. 

Following the standard procedure of QCD sum rules, we carry out the operator 
product expansion (OPE) of RelLjv^g 2 — ► —00) up to the lowest non-trivial order as 
is shown in Fig.l. We have neglected the perturbative a s corrections and took into 
account only the leading condensate (G 2 ). 



Fig.l 



To improve the OPE series, we apply the Borel transformation [14] to RelTy^g 2 ) 
which results in Ily y s{M 2 ) with M 2 being the Borel mass: 



2 ... it a s 2 , M 4 + 16m 2 M 2 — 16m 2 
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where m = ttT is the lowest Matsubara frequency, and 
9o(x) = 



ds^-e- s/M \ g x {x) = l~~ ds^~se- s/M \ (8) 



As for the imaginary part of Hs,v, we take 

lmU s (s) = a s 5(s-m s 2 ) + ^6(s-S s) (9) 

ImIV(s) = a v 5(s - m v 2 ) - — -=^- 9(s - S ov ), (10) 

4 vV 

where as (ay) is a constant and Sos (Sov) denotes the continuum threshold in the 
scalar (vector) channel. The structure of the continuum parts is simply obtained from 
Fig. 1(a) for large and positive q 2 . 

Eq.([|[7|) and (§,0) satisfy the Borel transformed dispersion relation 

B(M 2 ) = J lmU(s)e-^ds, (11) 
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The sum rule here are similar to that for the charm-quark system in (3+l)-dimensions. 
The charm-quark mass corresponds to our "Matsubara mass" (ttT) and the 4-dimensional 
vacuum condensate (GrL^o corresponds to our 3-dimensional condensate (Gjj). 
From (|IlD , one obtains the screening masses as 

(12) 



, m v 2 1 1 92i(p) - Mg^jg) - P 2 {^io(p) - ^io(P^)} + jfcggg! (13) 
27rT " ^io(p)-^io(^)- 1 t i ^ 5 (l + 2p 2 ) 

where p = 2m/M = 2nT/M, 

2 /"°° 2 

9ij(p) = 9i(p) -P 2 9j(p), g n (p) = 2e p x 2n e~ x dx, (14) 

and 

° " (2m) 4 " (2ttT) 4 ■ U5j 

Since the magnetic gluon condensate still exists above T c due to the non- 
perturbative effect in the space-like direct ion ||, we assume C ^ and introduce 
a scaling factor a relative to the vacuum magnetic-condensate in 4-dimensions, 

= h^) , (16) 

where (^B 2 ) = (1/2)(^G%,) - (1/2) (360 MeV) 4 0. Lattice QCD simulations 
suggests a ~ 1 for T c < T < (2 — 3)T C ||, while the free quark-gluon gas implies 
a = oo. 

Since the Borel mass is a fictitious parameter, the l.h.s. of eq . ([12|JT3D should be 
insensitive to the change of p in a certain window p min < p < p max (Borel window). 
We choose the window by the following conditions: (i) The first terms in the OPE 
(eqns.(|],[7D) are more than three times larger than the second terms so that the 
expansions are good enough, and (ii) The second terms in eqns.([],|II]) are more than 
three times smaller than the first terms so that final results are not sensitive to 
the detail structure of the continuum. The Sq's in eqns.(l9j,[r0f) are chosen so as to 



make (m sv /2irT) 2 least dependent on p within this window. The minima of m sv 
as a function of p, which we always find in the windows, are taken as the physical 
screening masses. Even if we adopt an average of msy{p) over the window as a 
physical screening mass, the results does not change much since the Borel curve is 
quite flat. 
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The gluon condensate and T 4 always enter as a ratio in C (see (fR]))- Therefore, 
as long as the good stability of m s v (p) is obtained as a function of p by the suitable 
choice of Sq's, an approximate scaling relation holds 

m sv {T ) a)^m sv {Ta 1 /\l). (17) 

Thus we can cover all values of a by rescaling T: The screening masses with large 
a (small magnetic condensate) at fixed T corresponds to that with a = 1 (standard 
value suggested in lattice simulations) at larger T. 

The final results for function of Ta ' are shown in Fig. 2. p depen- 

dence of m s at Ta 1 / 4 =0.4 GeV is also shown in Fig. 3 to show the good stability of 
our result in the Borel window. 



Fig.2 Fig.3 



From Fig.2, the following final results follow: 

(1) The screening masses in both channels deviate from the perturbative value 2ttT 

at moderate T and approach asymptotically to 2nT at large T. Furthermore, 
m v — m s ~ 300 MeV for wide range of T. These effects are induced by the 
non- vanishing magnetic condensate above T c . 

(2) The screening mass in the scalar channel is larger than 2ttT for the whole range of 

T. In fact, the coefficients of the gluon condensate C in the r.h.s. of eq.(|T2|) turn 
out to be positive in the Borel window, therefore the r.h.s. is always positive. 

The first feature (1), in particular, the positive gap m v — m s = 0(300 MeV) 
for wide range of T is quantitatively consistent with the result of the lattice QCD 
simulations [0. Although the lattice data are limited to relatively low T region 
(T < 2.5T C ) while our effective theory is valid at relatively high T > (2 — 3)T C ||, this 
agreement may give us a physical interpretation of the lattice data. The positive gap 
has been also obtained theoretically in ref . [[U], [TT[] where the spin-spin interaction in 



the "funny" space is essential. It is an open problem to understand the reason why 
our non-perturbative approach and their perturbative approach give similar result. 

The result (2) is qualitatively different from the lattice data JF| which show 
m s < 2ttT. Within the approximation adopted in this paper, we do not understand 
yet the origin of this difference. (Note that the funny space approach also shows 
m sv > 2nT if the spin-spin interaction is neglected |IT| .) So far, we have neglected 



several effects such as the a s corrections to Fig. 1(a) from the static modes (A,0) as 
well as the non-static modes. Since the effective quark mass m = nT is heavy at 
high T, the typical distance between q and q is small. Thus the above corrections, in 
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particular, the spin-independent attraction could give rise to non-negligible reduction 
of the screening masses just like the perturbative color-Coulomb interaction in the 
heavy-quark systems such as J/ip. QCD sum rules incorporating these effects are 
currently under investigation and will be reported elsewhere. 

In summary, we have studied effects of the non-perturbative magnetic conden- 
sate at high T on the hadronic correlations. The condensate gives rise to sizable 
corrections to the hadronic screening masses from their asymptotic value 2nT in a 
channel dependent way. The positive gap m v — m s and its magnitude are consistent 
with the current lattice data. 
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Figure captions 

Fig.l OPE diagrams in QCD3. (a) perturbative contribution and (b) the power 
correction by the magnetic gluon condensate. 

Fig. 2 The screening masses devided by 2nT against the scaled temperature Ta 1 / 4 
(GeV). 

Fig. 3 The behavior of {m s /2nTf for the scalar channel at Ta 1 / 4 = 0.4 GeV in the 
window of p defined in the text. 
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This figure "figl-l.png" is available in "png" format from: 
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